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Abstract 

We carried out a fluid dynamical simulation for a forced dripping faucet system using a 
new algorithm that was recently developed. The simulation shows that periodic external 
forcing induces transitions from chaotic to periodic motion and vice versa. We further 
constructed an improved mass-sprig model for the same system on the basis of data 
obtained from the fluid dynamical computations. A detailed analysis of this simple model 
demonstrated that the periodic motion is realized after a homoclinic bifurcation, although 
a stable periodic orbit is generated via a Hopf bifurcation which occurs just after a saddle 
node one. 

Chaos occurs widely in engineering and natural systems. In the past few years practical 
implementations of controlling chaos have been studied with great interest [1-6]. Such control 
schemes can be divided broadly into two categories, feedback control and nonfeedback control. 
The pioneering work on controlling chaos by Ott, Grebogi, and Yorke (OGY method) @] is a 
representative example of feedback control systems. The OGY method aims to stabilize one of 
many unstable periodic orbits embedded in the chaotic attractor through small time-dependent 
perturbations in an accessible system parameter. Although the OGY method is very general, it 
is difficult for some high-speed systems to implement the control procedure. In the nonfeedback 
control systems, on the other hand, the applied perturbation is independent of the system's 
state. So far, controlling chaos by applying a suitable weak periodic perturbation, which is 
sometimes called taming chaos [|J, has been studied in many chaotic dynamical systems (||. [5| 
and references therein). Recently, analyzing a constrained system in which a one-dimensional 
Poincare map is derived, Tamura et al. |4| showed that taming chaos occurs by a saddle 
node bifurcation. However, many issues in taming chaos are still not well understood, such as 
phase effect ||, phase diagram structures ||, and bifurcation structures for higher-dimensional 
dynamical systems f|. 
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In this Letter, we investigate taming chaos for a dripping faucet system. The dripping faucet 
without forcing has been studied intensively ([7-12] and references therein) since Shaw proposed 
a mass-spring model [[7] . Here we show from numerical simulations that periodic external forcing 
induces transitions between chaotic and periodic motion in this system. We further point out 
that a global homoclinic bifurcation is essential to realize the transition from chaotic to periodic 
motion. 

Before going to taming chaos, we present a counter phenomenon: chaos induced by a periodic 
perturbation, which was observed experimentally by Shoji fL2| |. By setting a periodically oscil- 
lating speaker on the top of the faucet, dripping time intervals {T n } between successive drops 
were measured. For sufficiently slow flow rates Q without any periodic perturbation the time 
series {T n } are almost equal at each flow rate (Fig. p] (a)). In contrast, the time series {T n } 
are distributed over a finite range for various flow rates when the speaker oscillates periodically 
(Fig. [I] (b)). Moreover, frequency entrainments are observed over a certain range of the flow 
rate. 

Similar results were obtained from our fluid dynamical computations (FDC) as shown in Fig. 
|2| (a)-(c). The algorithm of the FDC 0, which is based on Lagrangian description, was recently 
developed and succeeded in reproducing various complex behaviors observed experimentally in 
the dripping faucet under a constant flow velocity vq. In the present simulation, the flow 
velocity oscillates as 

v = v ^l + asm (^-^^j ■ (1) 
(The flow rate is defined as Q = ira 2 vo, where a is the faucet radius). Throughout the paper the 



units of the length, time, and mass are chosen as Iq = yT/pg (= 0.27 cm) , t = (T/pg 3 ) 1 ^ (= 
0.017 s) , and m = p/jj (= 0.020 g) , respectively, where T is the surface tension, p is the 
density, and g is the gravitational acceleration. (The numbers in parenthesis correspond to the 
water at 20 °C). The viscosity was rj = 0.002 (for the water) in units of 770 = (pT 3 /g) 1 ^. The 
FDC results (Fig. ^| (a)-(c)) qualitatively agree well with the experiment, although the number 
of the data points are not enough because of long computational time. 



Corresponding to the above FDC, we constructed an improved mass-spring model (see (10 
for details) with a periodic external force: 

d / dz\ dz . . /2ttA 

— \ m— \ = —kz — 'j- — hmo + Asin , (2) 

dt V dt J dt \ r J w 

dm 2 

— — = 7ra v = const., (3) 
dt w 

where z is the center of mass of the forming drop, and m is its mass. The damping parameter 
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7 was chosen as 7 = 0.008. The spring constant k depends on the mass as 

f 3.1m" - 58 (m < 0.584) 

k{m) = \ -31m 2 + 32m - 3.88 (0.584 < m < 0.891) . (4) 
[ (m > 0.891) 

Breakup of a drop is taken into account by assuming that a part of the mass is lost when the 
position z reaches a critical point z ciit . The remnant mass under the faucet is renewed as 

m r = Bm — C, when z = z cr i t , (5) 

depending on the total mass m at the breakup moment, where B = 0.068 and C = 0.053. The 
position and velocity just after a breakup moment are renewed as 



z = z 
z = 



when z = z crit . (6) 



The breakup position z crit and the renewed position zq were assumed as constant: z = 0.15 
and ;z cr it = 4. The relations (|J), (|| and the above parameters were obtained by referring to 
the results of the FDC for a = 0.183. 

When no external force is applied (A = 0), the remnant mass at the breakup moment, m x>n +ii 
is uniquely determined only by the previous remnant mass m r n , i.e., Poincare map is defined 
as m r>n _|_i = M(m In ), even though the return map of T n is not necessarily single- valued. In 



fact, This significant feature of the dripping faucet system was first found from the FDC |L0 
Under external forcing, on the other hand, the Poincare map just after the breakup moment 
is two dimensional because of an additional degree of freedom, the phase of the external force. 
The Poincare map is thus expressed as 

where 

9 = — modulo 1. 

T 

As shown in Fig. |2| (d)-(f), the mass-spring model can reproduce a bifurcation diagram in 
good agreement with the FDC (Fig. (a)-(c)). Moreover, the agreement with the experiment 
(Fig. [I] (a), (b)) is even better than for the FDC, because the numerical computations of this 
simple model can be carried out until transient behavior has completely vanished. 

Contrary to the above phenomenon, it is also possible to eliminate chaotic motion by applying 
a suitable periodic perturbation (taming chaos). An experiment by Katsuyama and Nagata for 
the faucet of about 5 mm in diameter, in which no periodic perturbation was applied, showed 
that a unit structure including chaotic transitions repeatedly appears in a bifurcation diagram 
in a wide range of the flow rate ||. A similar repeating structure in the bifurcation diagram 
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was obtained from both our FDC || and our improved mass-spring model [10]. Figure ^](a) is 



an example of a chaotic oscillation of T n obtained from the FDC for the faucet radius a = 0.916 
(2.5 mm) and the fluid velocity vq = 0.95, where no perturbation is applied. This chaotic 
oscillation is eliminated by a periodic perturbation (Eq. (j|D with a = 0.05 and r = 8.7) and a 
period-one (PI) motion with a period 2r is generated (Fig. 0(b)). 

In the mass-spring model, a PI motion (Fig. (|(b)) is also induced from a chaotic state (Fig. 
|4](a)) by periodic forcing. Referring to data obtained from the FDC for a = 0.916 (0.25 mm), 
parameters of the mass-spring model were chosen as (7 = 0.05, z cr i t = 5.5, z = 2.0, B = 0.2, 
and C = 0.3), and the spring constant was assumed as 

. , , f -11.4m + 52.5 (m<4.61) 

fc(m) = { (m>4.61) • (8) 

Figure § shows that, for the forcing amplitude A = 1.2 and the period r = 6, the system is 
entrained to the PI motion with the period 2r. 

We shall investigate how the entrainment is realized in the dripping faucet by analyzing the 
Poincare section of the mass-spring model in detail. Let x = (m r n , 6 n ) be a state point on the 
Poincare section just after the breakup moment (Eq. (0)), and let x = ^m r ,flj be a PI orbit 
satisfying x = M(x). If the PI orbit whose period is an integer multiple of the forcing period: 
T n = Nt = const., then the linear relation Eq. (|5|) yields an explicit expression for m r : 

_ BQNr + C 2 

m r = ^— — , Q = 7TCL Vq. (9) 

The PI orbits x should then satisfy 

xe/nM(/), (10) 

where I = {(m T ,6) \ m r = m r }, and the entire line I is mapped onto a curve M(J). For a 
specified value of the integer N, the PI orbits are obtained from Eq. (|9]). Figure |5](a) illustrates 
how the PI orbits with N = 2 are generated as the forcing amplitude A is increased, where / is 
shown by a dashed line and M(J) is shown by a solid curve. As A grows through A = 0.9347, 
/ intersects M(J). It was found numerically that two PI orbits xo and Xi, i.e., the intersection 
of / and M(I), are generated via a saddle node bifurcation where xi is a saddle. Despite the 
terminology 'saddle node', however, the other fixed point x is a source instead of a sink just 
after the saddle node bifurcation point. Therefore, there is only one attractor which is chaotic 
at this stage. It was found from a detailed analysis that bifurcations of the Poincare map of 
the present system can well be explained by a discretized version of a 2-dimensional flow which 
yields a homoclinic bifurcation: A scenario to a homoclinic bifurcation similar to that presented 
in Ref. |13| can be applied. Figure ^| illustrates invariant manifolds of the Poincare map (solid 
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curves with arrows) together with numerical examples of orbits near the manifolds (dots and 
thin lines). The saddle node bifurcation point is presented in Fig. ^(a). As the parameter A 
increases, the unstable fixed point xo is stabilized by a Hopf bifurcation (Fig. 0(b)). The basin 
boundary of xo is then an unstable limit cycle generated by the Hopf bifurcation. Further 
increase of A induces a homoclinic bifurcation (Fig. 0(c)). Up to this stage, the basin of x 
is a small closed region. After the homoclinic bifurcation, the basin of x , whose boundary is 
a heteroclinic orbit, is open. However, the chaotic attractor coexists and the motion remains 
chaotic because the basin of x is still small. After A passes through the critical value A = 0.963, 
the chaotic attractor disappears and the periodic motion xo is realized. The coexistence of the 
two attractors results in a hysteresis as shown in Fig. 0(a). If one starts from a large enough 
value of A (say A = 1.2, Fig. 0(b)) at which the PI motion is realized, and decreases the A 
value through A, the system keeps the PI motion till the Hopf bifurcation occurs at A ~ 0.936. 
At the critical point A = A, the system exhibits the following scaling property: 



n t ~ {A - A)- 1 - 1 for A > A. (11) 

Here n t is the average transient lifetime defined by the average iteration number before the 
orbit falls within a distance A from the stable PI orbit; the average was taken over many 
initial points on the Poincare section. The exponents and the value of A were determined by 
least square fitting. 

It is interesting that a similar structure of the Poincare section as in Fig. 0(d) was observed 



in a dripping faucet experiment by Pinto and Sartorelli JTTJ, although no periodic force was 
applied. The present analysis suggests that their result may possibly be explained if any 
oscillatory factor was generated inevitably in the experimental condition. 

In conclusion, we have demonstrated that periodic forcing induces transitions between chaotic 
and periodic motions in the realizable dripping faucet system. The dimension of the dynamical 
system generally increases by 1.0 when a forced perturbation is added. Then, a saddle node 
bifurcation that is not inherent in the original dynamics may occur. Higher-dimensional systems 
such as the dripping faucet may induce a global bifurcation further. It was shown that the 
change of the global structure in the phase space is essential in realizing the entrainment for 
the present system. Postnov et al. [TJJ] analyzed a periodically forced 4-dimensional flow, and 
found that a torus coexists with a stable periodic orbit between a saddle node bifurcation and 
a homoclinic one. In the present case, in contrast, chaos coexists with a stable periodic orbit 
between a Hopf bifurcation and a homoclinic one. 

We thank Professor M. Shoji for permission to copy his figure. This research was partly 
supported by Iketani Science and Technology Foundation. 
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Figure 1: Bifurcation diagrams in a dripping faucet experiment (a) without and (b) with 
periodic forcing. The faucet is 1.0 mm in inner diameter, 1.6 mm in outer diameter. In (b) the 
forcing period r = 0.05 s. [Reprinted by courtesy of Professor M. Shoji [12].] 
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Figure 2: Bifurcation diagrams (a)-(c) obtained from the fluid dynamical computations and 
(d)-(f) from the improved mass-spring model. (a)(d) the unperturbed system. The faucet 
radius a = 0.183 corresponds to 1.0 mm in inner diameter, r = 2 (= 0.034 s). 
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Figure 3: Time series of T n obtained from the fluid dynamical computations for (a) unperturbed 
and (b) periodically perturbed (a = 0.05, r = 8.7) system, a = 0.916, v = 0.95, N = 2. 




Figure 4: Time series of T n obtained from the improved mass-spring model for (a) unperturbed 
and (b) periodically perturbed (A = 1.2, r = 6) system, a = 0.916, v = 0.115, N = 2. 




Figure 5: (a) / is shown as a dashed line and their images M(J) are shown as solid curves at 
various forcing amplitudes A. (b) Orbits near period-one orbits on the Poincare section. 
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Figure 6: Invariant manifolds (solid curves with arrows) and examples of orbits (dots and thin 
lines). v = 0.115, r = 6. (a) A = 0.9347. (b) A = 0.940. (c) A = 0.9441. (d) A = 0.950. 
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Figure 7: (a) Log-log plot of the average transient lifetime n t vs (A — A), where A = 0.963. (b) 
Bifurcation diagram of T n with decreasing A. 
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